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1 Introduction 

The associated production of a Z/ 7 * or a 7 with one or more jets has been extensively studied in 
proton-proton collisions at the CERN LHC, by both the CMS [H-El and ATLAS ©IZl Collabora¬ 
tions. Precise measurements of these processes provide important tests of the Standard Model 
(SM) as well as crucial inputs in the determination of parton densities in the proton |PJ. Such 
measurements can improve the validation and tuning of the models used in Monte Carlo (MC) 
simulation. These processes are also important backgrounds in searches for new physics. 

In the limit of high transverse momentum of the vector boson V (pj) and at leading order 
(LO) in perturbative quantum chromodynamics (QCD), effects due to the mass of the Z boson 
(mz) are small, and the cross section ratio of Z -|- jets to 7 -|- jets as a function of pj is expected 
to become constant, reaching a plateau for pj > 300 GeV 0. (In this paper, production of 
Z/ 7 * -|- jets is denoted by Z -|- jets.) A QCD calculation at next-to-leading order (NLO) for 
pp Z -|-jets and pp ^ 7 -|-jets was provided by the BlackHat Collaboration ||9|. The NLO 
QCD corrections tend to lead to a decrease in the plateau value of the cross section ratio, while 
electroweak (EW) corrections are relatively small. However, at higher energies, EW corrections 
and QCD processes can introduce a dependence of the cross section on logarithmic terms of 
the form ln(pj/mz) that can become large and pose a challenge for perturbative calculations 
such as BlackHat+SHERPA llOl . A precise measurement of the (pp ^ Z-|-jets)/(pp —> 
7 -h jets) cross section ratio provides important information about the higher-order effects of 
these large logarithmic corrections at higher pj. 

In addition, searches for new particles involving final states characterized by the presence 
of large missing transverse energy (£t) and hard jets, as described for example in Refs. fTTlI 
and IT 2 I , use the 7 -|- jets process to model the invisible Z decays, Z —> vv, since the 7 -|- jets 
cross section is larger than the Z -|- jets process where the Z decays to leptons. Measurements 
of the cross section ratio for Z -|- jets and 7 -|- jets can help reduce uncertainties related to the 
Z —> pi7 background estimation in these searches. 

We present precise measurements of both production cross sections and the cross section ratio 
for these two processes as a function of pj. The results are compared with theoretical estima¬ 
tions. The data sample was collected at the LHC during the 2012 run with the CMS detector 
in proton-proton collisions at a center-of-mass energy = 8 TeV and corresponds to an inte¬ 
grated luminosity of 19.7 fb~^. The Z bosons are identified via their decays to pairs of electrons 
or muons. Measurements are made for different jet multiplicities (njets > 1/ 2, 3) and for a sub¬ 
set requiring a large hadronic transverse energy {Hj > 300 GeV, where Hj is the scalar sum 
of all selected jet pj after jet identification). These requirements are meant to mimic the phase 
space requirements for analyses searching for new physics with an all-hadronic signature. 

This paper is organized as follows: Section [^provides a brief description of the CMS detector; 
Section]^ gives details of the Monte Carlo generators used in this analysis; Section [^describes 
the event selection; Section|^contains details about the background subtraction and the unfold¬ 
ing of the detector effects; Section [^discusses the sources of systematic uncertainties; Section]^ 
presents the cross section and cross section ratio measurements for Z -|- jets and 7 -|- jets pro¬ 
duction. Section [^concludes with a summary of our results. 

2 CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal dia¬ 
meter, providing a magnetic field of 3.8 T. The CMS experiment uses a right-handed coordinate 
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3 Monte Carlo samples 


system, with the origin at the nominal interaction point, the x axis pointing to the center of the 
LHC ring, the y axis pointing up (perpendicular to the LHC ring), and the z axis along the 
anticlockwise-beam direction. The polar angle 6 is measured from the positive z axis and the 
pseudorapidity is defined as rj = — ln[tan(0/2)]. 

Within the superconducting solenoid volume are a silicon pixel and strip tracker, a lead tungstate 
crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter. 
Each subdetector is composed of a barrel and two endcap sections. Muons are measured in 
gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid. Extensive 
forward calorimetry complements the coverage provided by the barrel and endcap detectors. 
In the barrel section of the ECAL, an energy resolution of about 1% is achieved for unconverted 
or late-converting photons in the tens of GeV energy range. The remaining barrel photons have 
a resolution of about 1.3% up to a pseudorapidity of |?/| = 1, rising to about 2.5% at\r]\ = 1.4. 
In the endcaps, the resolution of unconverted or late-converting photons is about 2.5%, while 
the remaining endcap photons have a resolution between 3 and 4% |jl3|- The dielectron mass 
resolution for Z —)■ ee decays when both electrons are in the ECAL barrel is 1.8%, and is 2.7% 
when both electrons are in the endcaps. The electron momenta are estimated by combining en¬ 
ergy measurements in the ECAL with momentum measurements in the tracker [143 • Muons are 
measured in the pseudorapidity range \y\ < 2.4, with detection planes made using three tech¬ 
nologies: drift tubes, cathode strip chambers, and resistive-plate chambers. Matching muons 
to tracks measured in the silicon tracker results in a relative transverse momentum resolution 
for muons with 20 < px < 100 GeV of 1.3-2.0% in the barrel and better than 6% in the endcap. 
The px resolution in the barrel is better than 10% for muons with px up to 1 TeV IITSl . A more 
detailed description of the GMS system can be found in Ref. 111^ . 


3 Monte Carlo samples 

Monte Garlo simulation samples are used to correct the data for acceptance and efficiency and 
for the V -|-jets processes. They are also used to estimate the background to the Z -|-jets signal. 

The Z -|-jets signal is generated with the MadGraph (version 5.1.3.30) IfTTlI program. The 
leading-order multiparton matrix element (ME) calculation includes up to four partons (gluons 
and quarks) in the final state. The showering and hadronization of the partons, as well as the 
underlying event, are modeled by pythia (version 6.4.26) Ifl8l with the Z2* tune fl9 l. The kj 
MEM matching scheme 1201 with a matching parameter of 20 GeV is applied to avoid a double 
counting of final states arising in the ME calculation and the parton shower (PS). The events 
are generated with the GTEQ6L1 Il2ll parton distribution functions (PDE) and rescaled using 
a global next-to-next-to-leading-order (NNLO) R-factor to match the inclusive cross section 
calculated with FEWZ 3.1 |j22j. Backgrounds to Z -|-jets are generated using MadGraph with 
the same configuration as the signal events. These include top quark-antiquark pairs (tf) and 
EW backgrounds, such as W -|-jets and diboson processes (WZ, ZZ, WW). 

The Z -|-jets signal and tt background processes are also generated with SHERPA (version 1.4.2) 
[23], using the GTIO PDE Il24l . The cross section for the signal is also rescaled using a global 
next-to-next-to-leading-order (NNLO) fC-factor. Eor the tt background, the NNLO calculation 
provided by Ref. ||25l is used to calculate the NNLO R-factor. The background yields from EW 
processes are rescaled using the MGEM II26I NLO cross sections. 

In addition to these general purpose MG signal data sets, we use an NLO perturbative QGD 
calculation of Z -|- jets from the BlackHat Gollaboration S, which is available for a Z boson 
accompanied by up to three jets. These simulations use MSTW2008nlo68cl Il27] with ols = 0.119 
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as the PDF set, and the renormalization and factorization scales (jiR and Hf, respectively) are 
set to 

^R = fiF = H^ + Ej = EPt + \/?«|+ (pf)l 

i 

where pj is the transverse momentum of the ;th parton in the event and Hj is the scalar pj 
sum of all outgoing partons with pj > 20 GeV. The CTIO and NNPDF2.3 Il28l PDF sets with 
Ois = 0.119 are used as a cross check and to estimate the theoretical systematic uncertainties. 

The BlackHat+SHERPA simulated events are organized into different types of processes to 
facilitate the calculation. An NLO estimation at n jet level is obtained by combining tree-level 
(LO) calculations from the n + 1 jet case to n jet tree- and loop-level calculations. The Born 
and real emission calculations at both n and n + 1 jet levels are supplied by SHERPA, while 
BlackHat provides the NLO virtual loop-level correction terms. (For simplicity, Black- 
Hat+SHERPA will be referred to as BlackHat.) The structure of the generated files and the 
preselections used in the simulation, as well as more details about BlackHat, are described 
in Ref. [I0|. 

The 7 + jets signal is simulated by MadGraph, including up to four-parton final states in 
addition to the photon. Fixed-order cross section calculations for 7 -|- jets are affected by an 
instability due to dependencies on soft-gluon radiation, which can be overcome using all-order 
resummation |[29ll30ll . The background contribution due to multijets is determined using con¬ 
trol samples in data. The uncertainty in the photon purity is estimated with MG background 
samples simulated with PYTHIA6. Multijet events in the PYTHIA6 sample with signal-like be¬ 
havior can be enhanced by applying a filter that requires jet signatures with large electromag¬ 
netic deposits in the final state, e.g., jets with hadrons decaying into high-px photons. As an 
alternative method of estimating this background, we use a MadGraph sample that includes 
jet production with as few as two and as many as four outgoing partons in the ME calculation. 

We also simulate the 7 -|- jets signal using BlackHat. The overall procedure is analogous to 
the Z -|-jets BlackHat samples, and 7 -|- jets samples are available for 7 + 1, 2, and 3 jets. For 
7 -|- jets, we use the following renormalization and factorization scales: 

Pr — Pf — Efj -|- Ej = pj -|- pj. 

i 


The BlackHat production requires that the photons satisfy the Frixione cone isolation condi¬ 
tion II 3 TI 


Y,Ek@{S-Ri^)<n{5), 


for all 5 less than Sq around the axis of the photon. Here, Ri^ is the distance in p and azimuthal 
angle (f> between the f* parton and the photon, and 0 is the step function. The function T-L (S) 
is chosen such that it vanishes as ^ > 0. In particular. 


n {3) = E^e 


1 — cos 3 
1 — cos 3o 


n 


The Frixione cone in effect only adds contributions from partons which are within 3o of the 
photon. In the BlackHat samples, e = 0.025, 3o = 0.4, and n — 2. These were selected 
because a Frixione cone with these choices mimics the selections in the true on-shell photon 
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definition at particle level. Photon distributions using Frixione cone requirements are found to 
agree with those using cone isolation to within 1-2% 1321 . 

The simulation of the CMS detector is based on the Geant 4 package EH- The simulated 
events used for the detector level MC estimations are reconstructed following the same proce¬ 
dures used for the data. For our run, the average number of inelastic proton-proton collisions 
occurring per LHC bunch crossing was 21. The correct distribution of the number of pileup 
events overlapping the hard interaction process per bunch crossing is taken into account in 
the MC by reweighting the simulated minimum bias events to match the spectrum of pileup 
interactions observed in data. 


4 Event selection and object reconstruction 

The selection of Z -|-jets events begins by requiring two same-flavor high-px leptons (electrons 
or muons) at trigger level. The px threshold of the trigger objects is 17 GeV for the leading muon 
(the muon with the largest px) and 8 GeV for the subleading muon. The dielectron trigger re¬ 
quires the same thresholds of 17 and 8 GeV on the px of the leading and subleading electron 
candidates, respectively. Additionally, the trigger requires that the electron candidates be iso¬ 
lated from other energy deposits in the calorimeter so an isolation requirement is imposed on 
the electron track. 

Muons are reconstructed offline by a simultaneous fit of hits recorded in the silicon tracker 
and in the muon detectors HTSl . Electrons are reconstructed from energy clusters in the EGAL 
and tracking information IflTil . The two leading leptons are required to be of opposite electric 
charges and of the same flavor, with px > 20 GeV and \t]\ < 2.4. Eor both candidates, a match 
with a corresponding trigger object is required. The dilepton invariant mass, niu, is required 
to satisfy 71 GeV < niu < 111 GeV. This will be referred to as the "Z boson mass window". 

The particles in the event are reconstructed using the particle-flow (FE) technique Il34l l35l . 
which consists of identifying each single particle with an optimized combination of all subde¬ 
tector information. Depending on their signatures in the various subdetectors, particles fall into 
five different FE categories: muons, electrons, photons, neutral hadrons, and charged hadrons. 
The lepton candidates are required to be isolated from the other particles in the event, so to eval¬ 
uate the isolation a scalar px sum of FE objects is calculated in the cone AR = \/{Acp)^ -|- (Ap)^ 
around the direction of the object. The contribution from pileup to this isolation scalar px 
sum is subtracted using the average pileup energy per unit area in the rj-cp plane evaluated for 
each event Il36l . Eor electrons, the pileup-subtracted isolation sum is calculated in a cone of 
AR = 0.3 around the direction of the electron and is required to be below 15% of the electron 
Px. Eor muons, the radius is set to be AR = 0.4 and the isolation variable is required to be less 
than 12% of the muon px. Lepton reconstruction, identification, and isolation efficiencies are 
measured using the "tag-and-probe" technique as described in Ref. 1(371 . Efficiencies for sim¬ 
ulated events are corrected using p- and px-dependent scale factors to account for differences 
between data and simulation. Scale factors typically range between 0.98 and 1.02. 

The photons are reconstructed offline from energy clusters in EGAL IT^ . Events for the 7 -|- jets 
processes are selected at the trigger level, where the presence of a high-px photon candidate is 
required. Since the instantaneous luminosity of the LHG increased during the data-taking pe¬ 
riod, the threshold in px increased as well, and the lower px threshold triggers are prescaled 
in order to keep the rate at a reasonable level. An unprescaled trigger is available only for a 
photon with transverse momentum p^ > 207 GeV. In order to further reduce the rate, a loose 
shower shape cutoff < 0.24 is imposed at trigger level, where cr,jr; measures the extension 
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of the shower in pseudorapidity in terms of the energy-weighted spread within the 5x5 crystal 
matrix around the most energetic crystal in the photon cluster. For photon candidates, a match 
with a corresponding trigger object is required. For this analysis, only isolated high-pj photons 
located inside the barrel region of the detector (|p| < 1.4) are considered. We concentrate on 
photons inside the barrel region because the data size of the templates, described in Section 
allows for a precise purity determination in this region. Around 40% of the photons convert 
into e+e^ pairs inside the tracker material. Conversion track candidates are fitted from a combi¬ 
nation of ECAL seeded tracks and Gaussian sum filter Il38ll electron tracks originate from a com¬ 
mon vertex. The track pair is then matched to energy clusters in ECAL to identify a converted 
photon candidate. The final photon candidates are checked for possible overlap with electron 
candidates by looking for electron track seeds in the pixel detector or by using the character¬ 
istics of the track pair for converted photons. Isolation requirements are separately imposed 
on the pileup-corrected scalar pj sum of neutral and charged hadrons, as well as on additional 
photons inside a cone of AR = 0.3 around the photon candidate direction. MC over data scale 
factors for the selection efficiencies of unconverted and converted photons without the electron 
veto are measured using the tag-and-probe technique on Z —> e+e^ events; the scale factors for 
the electron veto efficiency on signal photons are determined using Z —> p+p “7 candidates. 
These scale factors range between 0.96 and 1.01 for photon candidates with pj > 40GeV. 

Jets are reconstructed from the four-momentum vectors of all PE objects. The anti-kx clustering 
algorithm Il39l is used here with a distance parameter of R = 0.5 in its FastJet BOl implemen¬ 
tation. The jets are clustered by four-momentum summation. The reconstructed PF candidates 
are calibrated separately to account for the nonlinear and nonuniform response of the GMS 
hadron calorimeter, especially for neutral hadrons. Gharged hadrons and photons are well 
measured in the silicon tracker and the EGAL, and therefore need only minimal corrections. 
Thus, the resulting jets require only small additional momentum adjustments. Jet energy cor¬ 
rections are obtained using Geant 4 simulated events generated with pythia6 . The energy 
contributions due to the presence of additional proton-proton interactions are subtracted from 
each jet using the measured pileup unit density in the event and the jet area 11361 . The p depen¬ 
dent corrections are adjusted using exclusive dijet events, while the pj-dependent corrections 
are adjusted using exclusive Z -|- 1-jet and 7 -|- 1-jet events in the data [41|. As a result of these 
adjustments, the reconstructed jets are corrected to the stable particle level ||42l. For PF jets, 
the jet energy correction factor typically ranges from 1 to 1.2. Jets originating from pileup are 
rejected using the criteria described in Ref. 1431 . This rejects 90-95% of pure pileup jets while 
keeping over 99% of jets from the primary interaction. Jet identification quality requirements 
are imposed in order to remove spurious jets caused by noise in the calorimeter. The remaining 
jets are accepted for the analysis if they satisfy pj > 30GeV and \p\ < 2.4. Additionally, jets 
within a radius of AR < 0.5 with respect to the axes of each lepton or photon candidate are 
removed. This cut affects a small number of jets. For both Z -|- jets and 7 -|- jets selections, the 
presence of at least one jet is required. 

The selection of Z -|- jets events is separate from the selection of 7 -|- jets events, and the two data 
samples are analyzed and corrected independently. The overlap between Z -|- jets and 7 -|- jets 
events is negligible. The analysis is repeated in four different, but not mutually exclusive, 
kinematic regions, with pj > 100 GeV and 

* ^jets ^ 1/ 

* ^jets ^ 2, 

* ^jets ^ 3, 

* Hi > 300 GeV. 


5.2 
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The rapidity of the Z boson is not restricted for the individual Z boson distributions. However, 
it is restricted to the rapidity range \y\ < 1.4 for the distributions of the ratio of pj to pj 
because the photon is measured only in this central rapidity range. Rapidity is defined as 
y = 2 In [(£ + pz) / (£ — Pz)]- The measured differential cross sections are binned in equal 
intervals of log^g pj{ GeV) of width 0.045 from 100 to 800 GeV, corresponding to the overlap 
region between the Z + jets and 7 + jets phase space. This binning ensures that as the number 
of events decreases, the bin width increases in a regular way. In terms of the photon purity 
determination (defined in Section [s!^ , the bins are chosen such that there are enough events in 
all bins in the final distribution to ensure a reliable measurement. 


5 Background determination and unfolding 

5.1 The Z + jets selection 

Events from the Z + jets process are selected as Z —> e+e^ and Z ^ candidates with one 

or more jets, as described in Section]^ The background-subtracted distributions are unfolded 
to the stable particle level for each decay channel separately and then combined. 

Several SM processes contribute to backgrounds to the Z -|- jets signal. For low p|, the most 
important background is tf production, whereas at higher pj values, diboson production is the 
dominant background. Gontributions due to W -|- jets and WW -|- jets are negligible for this 
analysis. The background contributions are subtracted using relative event rates predicted by 
MadGraph after an NNLO scaling for Drell-Yan and tt samples and an NLO scaling for the 
electro weak backgrounds. 

A cross check of the validity of the procedure for tt background estimation is performed using 
an ep control sample in data. This sample is largely dominated by tf production with an addi¬ 
tional contribution from fully leptonic Z —)■ decays. Both in absolute scale and shape, the 

simulation reproduces the dilepton transverse momentum spectrum (pj ) in the data within 
10%. This statement is valid both for a selection with a relaxed dilepton mass of > 60 GeV 
and a selection within the Z boson mass window as used for the final event selection. As a 
second check, the relative rate of ep events in data and MG are compared to those of dielectron 
or dimuon events as a function of the dilepton pj. Events with ep are selected by requiring the 
ep invariant mass to be either in the Z boson mass window or in the whole mass range. Events 
from e+e^ or p^p^ are selected in the Z boson mass window. All four distributions of these 
relative event rates from simulation are compatible with data within 10%. The tf background 
peaks at around pj ~ 100 GeV, where it amounts to 1.5% for the inclusive 1-jet selection and 
to 8 % for the inclusive 3-jet selection. In the high-Hx selection, it amounts to up to 12%. For 
Hj > 300 GeV, the relative rate drops below 0.5%, while in the ep channel no event is observed 
beyond p^ > 450 GeV. 

At around pj ~ 150 GeV, the EW background increases and reaches a plateau of about 5-7% 
for all phase space selections of the analysis beyond pj ~ 400 GeV. The rate of the com¬ 
bined EW backgrounds predicted by simulation in a control region of multilepton final states 
is checked with data in the following way. Instead of selecting the two leading leptons, and 
then enforcing the same flavor requirement, we instead select the first two leptons matching 
the trigger objects with the same flavor. The rate of events gained with respect to the baseline 
selection is largely dominated by diboson events. These additional data events are compared to 
estimations from MadGraph, finding an agreement within roughly 10% for all jet multiplicity 
phase space selections. This comparison is done in the range pj' < 300 GeV. 



5.2 The 7 + jets selection 
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The same selection criteria from data are used at the particle level: leading leptons are required 
to have px > 20GeV and |p| < 2.4, while jets are required to have pj > 30GeV within the 
region of |p| < 2.4. The particle level jets in simulation are obtained by clustering the gener¬ 
ated stable particles (after hadronization and including neutrinos) using the anti-Zcj algorithm 
with distance parameter of R = 0.5. Electrons and muons have different energy losses due to 
final state radiation at particle level. In order to compensate for these differences, we define 
a "dressed" level to make the electron and muon channels compatible to within 1%. This is 
achieved by defining in simulation a particle momentum vector by adding the momentum of 
the stable lepton and the momenta of all photons with a radius of AR = 0.1 around the stable 
lepton. All jets are required to be separated from each lepton by AR > 0.5. 

The background-subtracted detector-level distributions from data are unfolded to the particle 
level. The unfolding response matrix includes detector resolution effects and efficiencies. We 
use MadGraph to build a response matrix which allows us to map detector-level distributions 
to particle level. To quantify the bias introduced by the choice of the MG model, we use SHERPA 
as an alternative. The off-diagonal elements of the response matrices are small for both chan¬ 
nels. For the dielectron channel, 85-95% of all events in a given bin of the reconstructed p|® 
distribution are mapped onto the same bin at the particle level. For the dimuon charmel, at low 
Pj^ around 85% fall in this category, whereas at very high p^^, only 67% stay in the same bin 
at the particle level. The remaining events typically fluctuate to directly neighboring bins. The 
iterative method used by d'Agostini 1441 , as implemented in the RooUnfold package Il45l , is 
used to regularize the inversion of the matrix. Subsequently, the unfolded distributions from 
Z —)■ e+e“ and Z —)■ are found to be compatible. They are combined using the best 

linear unbiased estimator Il46l to obtain the final distributions. The resulting averaged leptonic 
Z -h jets distributions from both charmels are not corrected to the total cross section. 

5.2 The 7 + jets selection 

After selecting 7 -|-jets events (as described in Section]^, the photon signal purity is determined 
in each px bin. The main background is due to QGD multijet production, where either one of the 
jets, or an electron or tt® inside a jet, is misidentified as a photon candidate. Since simulations 
do not provide a reliable description of this background, the purity, which is defined as the 
number of true isolated photons from the hard scattering versus the number of all photon 
candidates, is determined from data. At the particle level, a true isolated photon is defined 
as a prompt photon, around which the scalar sum of the px of all stable particles in a cone 
of radius AR = 0.4 is less than 5 GeV. Similarly, at the detector level, for each px bin of the 
photon spectrum, the purity is determined through a fit of the photon isolation sum variable 
^ph' defined as the scalar px sum of all other PF photons around the axis of the selected photon 
candidate, inside a cone of AR = 0.4. The sum is corrected for the pileup contribution and the 
energy deposit ("footprint") of the selected photon candidate itself. The distribution for the 
data is fitted as a sum of signal and background template distributions, in each px bin, in order 
to calculate the purity/: fp^(data) = / Jp^(signal) + (!—/) fp^(background). 

In order to model the contribution of the underlying event to the photon component isolation 
sum around the signal photon candidate, a signal template is obtained from the data through 
the random-cone (RG) method H47l . After selecting photon candidates fulfilling a requirement 
on the shower shape of (7^,^ < 0.011, a cone of AR = 0.4 is randomly chosen in (p at the same 
rj as the photon candidate, excluding the back-to-back direction to avoid selecting any recoil¬ 
ing jet. The candidate cone is rejected if it contains objects originating from a hard interaction, 
e.g., jets with px > 30 GeV or photons with px > 20 GeV. The RG templates show a very 
good agreement over orders of magnitude between data, simulation, and the true-photon MG 
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templates obtained by matching a detector-level photon candidate with an isolated photon at 
the particle level. Background templates are constructed by selecting photon candidates in the 
data with an inverted shower shape requirement, 0.011 < CTrjri < 0.014. Since there are a small 
number of background events with high pj photon candidates, the templates are obtained in 
wider bins of Pj than used in the final analysis. After construction of the templates, a binned 
maximum-likelihood fit to the Jp^ data distributions is performed as a sum of the signal and 
background template distributions. The statistical uncertainty on the fit includes the effect of 
the limited template sample size. An example of the fit for the photon component of the pho¬ 
ton isolation can be seen in Fig. [^(left). The results of such fits are displayed in Fig. (right), 
where the black dots represent the measured purity and the solid lines represent the statistical 
uncertainties including the template uncertainties. The purity ranges from around 65% in the 
lower Pj bins to 90% at high pj. The shaded band represents the total statistical and systematic 
uncertainties in the purity measurement. The difference in the fractions obtained from the clo¬ 
sure test on simulation is treated as the systematic uncertainty. True MC signal (background) 
templates are determined using identified photons at detector level matched (not matched) to 
a particle-level photon. The closure test uses the same approach for deriving templates that 
was used in data: it takes simulated samples and compares the resulting templates of the data- 
driven approach with the templates constructed exploiting the MC-truth information. This is 
the dominant source of systematic uncertainty in the purity estimate. Other effects, such as a 
change in the (Ttjtj requirement, are found to be negligible. The systematic uncertainties of the 
purity estimate are discussed further in Section]^ The data yields in each pj bin, after correct¬ 
ing for purity, are unfolded to the particle level with a procedure identical to the one used for 
the Z -|-jets process. Over the whole pj spectrum, the diagonal elements of the response matrix 
contain more than 90% of the events. 



Figure 1: The purity fit on the photon component of the photon isolation in the photon 
transverse momentum bin between 100 and 111 GeV in data (left). The photon purity as func¬ 
tion of the photon transverse momentum (right). The dots are the data points, the dot-dashed 
line is the signal template, and the dotted line represents the background component. The solid 
line represents the fit and the legend shows the resulting purity fraction of 66.7%. 
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6 Systematic uncertainties 

For Z + jets and 7 + jets we consider the following uncertainties: the jet energy scale (JES) un¬ 
certainty, the jet energy resolution (JER) uncertainty, uncertainties due to the MC model from 
the unfolding procedure (UF), the pileup uncertainty (PU), and the luminosity uncertainty 
(Lumi). Systematic effects specific to 7 -|-jets are those related to the photon energy scale (7 
ES) and the purity determination (7 Pur). For Z -|-jets, we consider the background subtraction 
(BG) and the lepton (muon and electron) momentum scale (LS) and resolution (LRES) uncer¬ 
tainties, as well as uncertainties in the lepton efficiency and isolation. We also consider lepton 
efficiency scale factors (lep SFs) for the Z -|- jets events. For both processes, the uncertainty 
associated with the luminosity is a flat 2 . 6 % over the whole range P48j. 

The systematic effects due to scale uncertainties affect the data, so we vary the momenta of 
the jets or the leptons independently within their uncertainties and rerun the unfolding on the 
shifted distribution. The differences in the final results are taken to be the systematic uncer¬ 
tainties. 

The uncertainty due to the JES affects distributions through the jet pj threshold Ell. For central 
jets, the JES uncertainty is around 3% at 30 GeV, decreasing to 1% at 100 GeV. Therefore, in the 
'^jets = 1 case, the boson pj is almost completely unaffected. For the 2- or 3-jet inclusive phase 
space selection, the requirement of additional jet activity increases the uncertainty in the photon 
and Z spectra due to JES to 5-10% over the whole pj range. In the Hj > 300 GeV selection, the 
JES uncertainty is around 5-7% at low pj and below 1% for pj > 400 GeV. 

Systematic effects due to purity and background subtraction have to be applied prior to un¬ 
folding to evaluate the uncertainty. The other uncertainties affect the response matrix in the 
simulation, and the unfolding is performed with these modified matrices to determine the rel¬ 
ative uncertainty. For example, we modify the resolution of jets in the MG and then calculate 
a new response matrix with these modified resolutions. The difference between this result and 
the nominal result is taken as the uncertainty due to JER. 

Eor the Z -|-jets process, the dominant sources of uncertainty are the lepton SEs and the LS in the 
^jets > 1 case and the JES uncertainty otherwise. The uncertainty due to the background sub¬ 
traction is typically below 1%. The lepton resolution uncertainty has an effect that is typically 
less than 0.5%. The effect of the electron energy scale uncertainty increases with pj from 1% 
at 40 GeV to 5% at 800 GeV. Eor muons, the scale uncertainty has an effect <1% up to 250 GeV, 
which increases up to 15% at high pj. Above 200 GeV, the track becomes very straight, and 
so the influence of the muon system becomes more relevant with respect to the tracker for the 
muon Pj distribution. This leads to an increase in the muon scale uncertainty. 

Eor the unfolding procedure, an additional check using the matrix obtained from SHERPA in¬ 
stead of that from MadGraph is performed, resulting in a cross section uncertainty of 2 - 3 % for 
all phase space regions. A cross check using the Singular Value Decomposition regularization 
method t 49l for the unfolding shows negligible deviations. 

The JER is measured to be about 5% larger than predicted in simulation for the central detector 
part < 1.4) with an uncertainty of about 5%, and roughly 10% larger in the endcaps with 
an uncertainty of roughly 7% BTl . The JER and pileup uncertainties in Z -|-jets events have 
values typically below 1%. The uncertainties on Z -|- jets are summarized in Table 

Aside from the pj spectrum, we also consider two additional variables for the Z -|- jets final 
state: the ratio of the pj to Hj and to the log^Q / Px)' where pj is the transverse momen¬ 
tum of the largest jet in the event. Most uncertainties in these distributions are similar to those 
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7 Results 


described above for the pj spectra, with the exception of the JES uncertainty. The latter has a 
larger influence on hadronic quantities {Hj, pj of the jets, and njets) which enter the distribu¬ 
tions directly, rather than through phase space selections. 


Table 1: Systematic Uncertainties for the pf Spectrum 


Process 

JES 

(^jets ^1) 

JES 

(otherwise) 

JER Lep SFs 

UF 

PU, BG 
LRES 

LS 

Lumi 

Z —> e+e^ 

1-3% 

5-10% 

<1% 3-4% 

2-3% 

<1% 

1-5% 

2.6% 

Z — 

1-3% 

5-10% 

<1% 2.5-5.5% 

2-3% 

<1% 

<1% 

2.6% 


For the 7 + jets process, the dominant uncertainty is due to the photon purity. This is a result 
of the difference between the shapes of the templates defined using the data-driven techniques 
from above and the distributions of the "true" templates for isolated photons in simulated 
events. Data samples are generated using the distributions of isolation variables for every bin 
of each variable with the fractions measured in data. Each of these is fitted with templates built 
in MC using the same techniques as on data, and the average difference between these fitted 
fractions and the generated fractions is quoted as the systematic uncertainty due to photon 
purity estimation Il47ll . This difference is around 10% when pj « 100 GeV and it decreases to 
roughly 4% at pj ~ 400 GeV for the inclusive njets > 1 selection. A change in the selection 
criteria on (7,jrj leads to negligible effects on the purity estimation. The background templates 
do not show any dependence on Hj or on the number of jets in the analysis. Therefore, the 
same background templates of the inclusive selection are used for all phase space regions and 
a similar template uncertainty is obtained. 

The JER uncertainty has a negligible effect in the analysis region pj > 100 GeV. For the 2- or 
3-jet phase space and high Hj > 300 GeV selection, the resolution uncertainty has an effect 
around 0.5-1.5%. The effect of the 7 ES uncertainty on the cross section measurement is con¬ 
stant across the whole range and less than 3%. The unfolding uncertainty is estimated by using 
an unfolding matrix from PYTHIA6 simulation and is around 2%. The uncertainty in the pileup 
interactions is evaluated by rescaling the cross section of minimum bias events by 5% in the 
MG reweighting procedure. Typically, these uncertainties are very small, below 0.5%. These 
systematic uncertainties for the 7 -|-jets process are summarized in Table 

Table 2: Systematic Uncertainties for the pj Spectrum 


Process 

JES 

(^jets ^1) 

JES 

(otherwise) 

JER UF 

PU 7 Pur 

7 ES Lumi 

7 

1-3% 

5-10% 

0.5-1.5% 2% 

<0.5% 4-10% 

3% 2.6% 


7 Results 

7.1 Differential cross sections 

In Figs. 0 and we present the measured differential cross sections as functions of the pj 
and the p^ for two selections of Z -|- jets and 7 -|- jets events (njets > 1 and njets > 2) and com¬ 
pare them with estimates from BlackHat and MadGraph+PYTHIA6 . In Fig. we present 
the ratio of the inclusive 2-jet events to the inclusive 1-jet events. For Z -|- jets, we also com¬ 
pare the data to SHERPA results. The NFO BlackHat estimate is corrected for nonpertur- 
bative effects (hadronization and underlying event) using MadGraph+pythia6 . These cor¬ 
rections are typically around 2%. We use the n-jets BlackHat sample for comparison with 
data and other MG generators in the corresponding inclusive n-jets selection. The Z -|-jets sim¬ 
ulations from SHERPA and MadGraph+pythia6 are rescaled by a constant NNFO R-factor 
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of K = 1.197, as calculated with FEWZ 3.1 Il22l . while for 7 + jets the LO cross section from 
MadGraph+pythia6 is used as no NNLO fC-factor is available for 7 + jets. In all figures, 
the hatched band surrounding the data points represents the total uncertainty in the measure¬ 
ment, while the error bars show the statistical uncertainty. Similarly, in the MC/ data ratio plots, 
the error bars around the points centered at one represent the relative statistical uncertainties 
on the data, while the hatched band represents the relative total uncertainty of statistics and 
systematics on the data. The shaded bands around the MC simulation/data ratios for Mad- 
Graph+pythia6 and SHERPA represent the statistical uncertainty (stat. unc.) in the simulation. 
The outer hatched band around the BlackHat/ data ratio (using MSTW2008) shows the total 
uncertainty of the estimate due to PDF and scale variations, while the inner hatched band in¬ 
dicates the uncertainty due to the variations within the MSTW2008 eigenvector set KSOl . Anal¬ 
ogous variations using the CTIO and NNPDF2.3 PDF sets lead to similar uncertainties. Not 
shown in the figures is the statistical uncertainty for the BlackHat calculations that amounts 
to less than 1-3% for njets >1/2 and to 5-10% for njets > 3 in the pj spectra. In the distribu¬ 
tions of the observables pj/Hj and log^Q (j’j/pj^, the statistical uncertainty is 6% except in 

the tails where there are fewer events. In the pj spectra, the statistical uncertainty is 3-5% in the 
^jets > 1 case and 4-10% in the njets > 2 and 3 cases. The fluctuations seen in the BlackHat 
distributions between adjacent bins are statistical in nature. Overlaid are BlackHat estimates 
using the NNPDF (dashed) and GTIO (dotted) PDF sets. 

In the Z -|- jets distributions for both phase space selections {njets > 1 and njets > 2, Figs, [^an d|3), 
we observe the same qualitative behavior of the ratio of the MadGraph+pythia6 simulation 
to data, which is flat about unity up to around 150-200 GeV and then increases to about 1.3 at 
higher pj. Estimates from SHERPA are lower than the data for pj < 50 GeV, while for higher 
Pj they increase to around 20% higher. In the njets > 1 case, BlackHat shows a flat ratio with 
respect to data starting around pj ~ 100 GeV, but underestimates the yield seen in data by 
8-10%, whereas in the njets > 2 case, BlackHat agrees with the data within the uncertainties 
for the whole range. For all multiplicity phase space selections, the systematic uncertainty in 
the MSTW2008 PDF set is 2-3% in the BlackHat estimate. The central points of GTIO show 
a difference compared to MSTW2008 of at most 4%, whereas NNPDF shows a variation of 
2%. The scale uncertainty for MSTW2008 in the BlackHat estimate, as determined through 
independent variations of the renormalization and factorization scales by factors of 2 and 0.5, 
leads to an envelope with values of typically 5-10%. 

In the njets > 2 and njets > 1 7 -|-jets case, BlackHat reproduces the shape of the data distri¬ 
bution, but underestimates the rate by approximately 10 - 15 % throughout most of the range. 

In Fig. 1^ we see that the inclusive 1-jet over inclusive 2-jet pj cross section ratio increases 
until a plateau is reached at around 350 GeV. The systematic uncertainties are treated as fully 
correlated in the ratio. The distributions are well predicted by MadGraph+pythia6 in both 
channels. SHERPA underestimates the relative rate of inclusive 2-jets events. For BlackHat, 
the inclusive 2-jet generated sample is used to predict the 2-jet rate and to compute the ratio 
with the predicted rates from the inclusive 1-jet sample; BlackHat overestimates the ratio by 
10% for pj > 100 GeV in both the 7 -|- jets and Z -|- jets cases. 

For Z -|-jets, we study the variables pj/Hj, shown in Fig. and log^g (pj/Pj'jf shown in 
Fig- HI which allow us to test the validity of NFO estimations. In particular, we examine these 
distributions as quantities where NFO estimations might reach their calculational limit due to 
large logarithms or where missing higher-order effects could play a larger role. 

For events which contain a dominant high-px jet, Pj/Hj tends to unity as the jet carries most 
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19.7 fb ' (8TeV) 




p \ [GeV] 


Figure 2: Top left: Differential cross section for Z boson production as a function of pj for 
an inclusive Z + jets, njets > 1 selection of detector-corrected data in comparison with esti¬ 
mations from MadGraph+PYTHIA6 , SHERPA, and BlackHat. Top right: Differential cross 
section for photon production as a function of pj for an inclusive 7 -|- jets, njets > 1 selection 
for central rapidities |y'^| < 1.4 in detector-corrected data is compared with estimations from 
MadGraph+pythia6 and BlackHat. A detailed explanation is given in Section 7.1 The 


bottom plots give the ratio of the various theoretical estimations to the data in the Z -|- jets case 
(bottom left) and 7 -|- jets case (bottom right). 
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Figure 3: Top left: Differential cross section for Z boson production as a function of pj for 
an inclusive Z + jets, njets > 2 selection of detector-corrected data in comparison with esti¬ 
mations from MadGraph+PYTHIA6 , SHERPA, and BlackHat. Top right: Differential cross 
section for photon production as a function of pj for an inclusive 7 -|- jets, njets > 2 selection 
for central rapidities |y'^| < 1.4 in detector-corrected data is compared with estimations from 
MadGraph+pythia6 and BlackHat. A detailed explanation is given in Section 7.1 The 


bottom plots give the ratio of the various theoretical estimations to the data in the Z -|- jets case 
(bottom left) and 7 -|- jets case (bottom right). 
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Figure 4: Ratio of the inclusive rates for njets > 2 and njets > 1 versus the transverse mo¬ 
mentum of the boson for Z -|- jets in detector-corrected data compared to estimations from 
MadGraph+pythia 6 , SHERPA, and BlackHat (top left) and for 7 -|-jets for central rapidities 
ly'^l < 1.4 in detector-corrected data compared with estimations from MadGraph+pythia 6 


and BlackHat (top right). A detailed explanation is given in Section 7.1 The bottom plots 


give the ratio of the various theoretical estimations to the data in the Z -|- jets case (bottom left) 
and 7 -|-jets case (bottom right). 
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of the Pt of the event. Events that populate the high-end tail of the distribution have either 
additional jets outside of the acceptance in the forward region or additional hadronic radiation 
that is not clustered in jets with > 30 GeV. In hadronic searches for new physics, these 
events contribute to signatures with a high Ex /Hj ratio. Almost all events with two or more 
jets inside the jet acceptance selections have pj/Hj values below one. This behavior can be 
observed in Fig. increasing the number of required jets leads to a shift of the complete dis¬ 
tribution towards lower values. The nonperturbative corrections are slightly larger, typically 
below 5% in the bulk of the distribution, reaching 10% in the tails for all variables examined 
here. Overall, MadGraph+pythia6 predicts the rate and shape best up to the tails, while 
SHERPA shows differences in both shape and rates. The BlackHat generator performs well 
for the bulk of the distribution, but fails to reproduce the tails. This is especially evident in the 
high-end tail of the distribution, where we see a sharp drop in the ratio of BlackHat simula¬ 
tions to data. In this portion of phase space, BlackHat is effectively reduced from an NLO to 
LO calculation as the n+1 jet LO calculation in the inclusive n jet case dominates here, whereas 
the other portions provide negligible contributions. This feature is also confirmed by the sharp 
increase of scale uncertainties in BlackHat estimates, which have a step-like increase from 
below 10% to around 60% at this point. Therefore, this sharp change in the BlackHat over 
data ratio (e.g., around 1.2 and 1.1 in Fig.|^ is expected and indicates the "boundary" between 
the regions where a fixed-order calculation gives a suitable estimation and where we would 
need the parton showering to add soft jets or jets in the forward regions of the detector. Ad¬ 
ditionally, we can use this to check for any large logarithmic contributions in the lower end of 
this pj/ Hj distribution. We see from the agreement in both the 2- and 3-jet cases that there is 
no evidence of any such contributions. 


The distribution of the second variable, log^Q {pj/pj), shown in Fig. 


shows similar behavior. 


For events with exactly one jet, the Z boson and the jet are back-to-back, with pj ^ pj, and the 
distribution peaks around zero. Events where the Z boson is the dominating object will have 
positive values. If the Z boson carries less pj than most of the jets, the variable has negative 
values. With increasing jet multiplicity the distribution still peaks around zero, but broadens. 
Figure]^ shows a comparison of MG estimates to data, which is unfolded to particle level. The 
MadGraph+pythia6 calculation performs well in estimating the behavior of the data for all 
inclusive multiplicty selections, but there is a slope within uncertainties in the MG/data plot. 
On the other hand, BlackHat performs well in the middle range, but the behavior in the tails 
indicates that jet production due to higher-order diagrams is missing. In the Z -|- jets, njets > 2 
phase space we observe a drop at the value log^g(pj/ pj) = 0.3 ~ log^Q(2), corresponding to 
pij! « with both jets recoiling against the Z boson direction. The distribution drops at the 
point where the third-leading jet becomes relevant. Since we use the inclusive 2-jet Black¬ 
Hat sample in that phase space, 3-jet events are only available as FO contributions in the real 
part. Therefore, the estimation is effectively an FO calculation at that point onwards, and sub¬ 
sequently becomes less precise and the scale variation uncertainty increases to around 30% at 
that point. 


7.2 The Z /7 ratio 

In order to compare the cross sections for Z -|- jets and 7 -|- jets, the rapidity range of the bosons 
is restricted to \y^\ < 1.4 because this is the selected kinematic region for the photons. The 
ratio of the differential cross sections as a function of px is measured in the four phase space 
regions: njets > 1/ 2, 3, and Hj > 300 GeV, njets > 1- 

Statistical uncertainties in the ratio are propagated using the diagonal terms of the covariance 
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Figure 5: The measured distribution of the observable pj/Hj ratio for njets > 2 (top left) 
and njets > 3 (top right) for Z + jets in detector-corrected data compared with estimations 
from MadGraph+pythia6 , SHERPA, and BlackHat. A detailed explanation is given in Sec¬ 
tion 7.1 The bottom plots give the ratio of the various theoretical estimations to the data in the 
^jets > 2 case (bottom left) and njets > 3 case (bottom right). 
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Figure 6: The measured distribution of the observable logjg pj/pj ratio for njets > 2 (top left) 
and njets > 3 (top right) for Z + jets in detector-corrected data compared with estimations 
from MadGraph-i-pythia6 , SHERPA, and BlackHat. a detailed explanation is given in Sec¬ 
tion 7.1 The bottom plots give the ratio of the various theoretical estimations to the data in the 
'^jets > 2 case (bottom left) and njets > 3 case (bottom right). 
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matrices. The sources of systematic uncertainty such as the JES, the luminosity uncertainty, and 
the JER are correlated between Z + jets and 7 + jets and therefore cancel in the ratio. The re¬ 
maining uncertainties are results of the photon purity measurement, unfolding uncertainty, the 
uncertainties in the efficiency determination for photons, and the lepton energy or momentum 
scale uncertainty. 

The resulting ratio distributions are shown in Eigs.[^and|^for all selections. The Z -|-jets selec¬ 
tion with the requirement Hj > 300 GeV enhances the presence of events with large hadronic 
activity. 

In all phase space regions, we observe a ratio which saturates around pj ~ 300 — 350 GeV. This 
agrees with the LO estimations stating that the main distinction between the two processes is 
the mass difference, with the second difference being the different couplings. 

In the inclusive njets > 1 selection (Eig.j^, the plateau value is 


o'z^£+e-{PT > 314 GeV) 
dyipj > 314 GeV) 


= 0.0322 ± 0.0008 (stat) ± 0.0020 (syst). 


( 1 ) 


Here Rauep is the plateau value of the ratio of the dilepton Z cross section and the 7 -|- jets cross 
section for the last seven bins {pj > 314 GeV). This translates into the ratio of the total cross 
sections of Rtot = 0.957 ± 0.066 when divided by the average leptonic branching fraction of 
(3.3658 ± 0.0023)% fSH. 

The estimation from MadGraph+PYTHIA6 is overlaid in Eigs. and where the LO esti¬ 
mation is used to compare Z -|- jets and 7 -|- jets differential cross sections at the same order of 
perturbative expansion. Although MadGraph+pythia6 does not reproduce the high-px tail 
for either Z -|- jets or 7 -|- jets, the shapes of the curves are similar for both processes and their 
ratio is flat. Using LO cross sections, MadGraph+pythia6 predicts a ratio with a value of 
Rmg = 0.0391, which is higher than that observed in data by a factor of 1.21 ± 0.08 (stat+syst). 
No clear trend away from a flat ratio is observed. Higher-order effects beyond LO, which could 
lead to a rise or fall in the plateau region, are smaller than the experimental uncertainties. 

The BlackHat estimation is also overlaid in Eigs. |^-[^and reproduces the 1-jet and 2-jet ratio 
to within 10% across the entire range. It reproduces the Hj > 300 GeV case accurately in the 
low-px regime and results in an approximately 20% overestimation in the high-px range. In 
the region where pj < 300 GeV, the scale uncertainty grows to roughly 30%. This corresponds 
to the region where BlackHat fails to reproduce the pf and pj spectra separately. Inclusive 
fixed-order calculations are not designed to model this selection of high jet activity with a com¬ 
paratively low boson px. In the 3-jet case, BlackHat overestimates the ratio by approximately 
25%, but agrees with data starting around the plateau region of approximately 300 GeV. 

We calculate the scale and PDE uncertainty bands for BlackHat using the scale and PDE 
uncertainty envelopes from the pj and pj spectra. If we correlate the different renormali z ation 
and factorization scales (pR and pp), the envelope decreases to approximately 2 %, whereas if 
we take the scales as completely anticorrelated, we see a band of approximately 10 % in the 
bulk. However, we know that the former underestimates the theoretical uncertainty due to 
renormalization and factorization scales, and the latter overestimates it. The estimation of this 
uncertainty has been discussed in the literature, and has been examined by comparing different 
theoretical computational estimations (|9j and Il32ll l. Both of the previously mentioned methods 
misrepresent the actual uncertainty due to the renormalization and factorization scales. We 
therefore choose the larger relative scale uncertainty band from each process as an estimate of 
the uncertainty on the final ratio. Using the NLO cross sections, BlackHat predicts the Rdiiep 
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ratio with a value of Rbh = 0.03794, which is higher than that observed in data by a factor of 
1.18 it 0.14 (stat + syst). 



19.7 fb ' (8TeV) 


1.4 

CD 

§1.2 

5 1.0 

« 0.8 
CQ 

0.6 
CD 1-4 

CD 

Q 1.2 

2 1-0 
o 

^ 0.6 



4 Scale 



I MadGraph stat. unc. 


100 200 300 400 500 600 700 


800 


p;'lGeV] 



100 200 300 400 500 600 700 800 


pf[GeV] 


Figure 7: Differential cross section ratio of averaged Z —)■ (e^e^ + }i^}i~) over 7 as a function 
of the total transverse-momentum cross section and for central bosons < 1.4) at different 
kinematic selections in detector-corrected data. Top left: inclusive (njets > 1); top right: Hj > 
300 GeV, njets > 1. The black error bars reflect the statistical uncertainty in the ratio, the hatched 
(gray) band represents the total uncertainty in the measurement. The shaded band around the 
MadGraph+pythia6 simulation to data ratio represents the statistical uncertainty in the MG 
estimation. The bottom plots give the ratio of the various theoretical estimations to the data in 
the njets > 1 case (bottom left) and Hj > 300 GeV case (bottom right). 


8 Summary 

Differential cross sections have been measured for Z -|- jets (with Z —>• £^£^) and isolated 
7 -h jets as a function of the boson transverse momentum, using data collected by GMS at 
= 8 TeV corresponding to an integrated luminosity of 19.7fb^^. The estimations from the 
MG multiparton LO+PS generators MadGraph+pythia6 and SHERPA have been compared 
to the data. We find that the pj spectra for Z -|- jets and 7 -|- jets are not well reproduced by 
these MG models. We observe a monotonic increase of the MG simulation/data ratio with in¬ 
creasing vector boson pj. Using the NLO generator BlackHat simulation, we find a smaller 
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Figure 8 : Differential cross section ratio of Z —> (e+e^ + ) over 7 as a function of the total 

transverse-momentum cross section and for central bosons (|i/^| < 1.4) at different kinematic 
selections in detector-corrected data. Top left: 2-jet (njets > 2); top right: 3-jet (njets > 3). The 
black error bars reflect the statistical uncertainty in the ratio, the hatched (gray) band represents 
the total uncertainty in the measurement. The shaded band around the MadGraph+pythia6 
simulation to data ratio represents the statistical uncertainty in the MC estimation. The bottom 
plots give the ratio of the various theoretical estimations to the data in the njets > 2 case (bottom 
left) and njets > 3 case (bottom right). 
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discrepancy in shape between data and simulation, indicating that it is likely related to missing 
higher-order effects. 

We have also studied the distribution of the ratios of pf and hadronic quantities (Hj and pj) 
in Z -|- jets. We find that these agree with the LO+PS estimation over the whole range when an 
NNLO K-factor is applied. The NLO BlackHat estimation is accurate in a subrange where 
the NLO estimation is expected to perform well. 

In addition, we presented a measurement of the ratio of the Z + jets to 7 + jets cross sections in 
four phase space regions: n^ets > 1/ 2, 3, and Hj > 300 GeV, njets > 1- MadGraph+pythia6 
(LO+PS) overestimates the data by a factor 1.21 ± 0.08 (stat+syst), whereas BlackHat (NLO) 
overestimates the data by a factor 1.18 ± 0.14 (stat+syst) in the plateau region i.e., for pj above 
approximately 300 GeV. As a function of the vector boson transverse momentum, these factors 
are at similar values of around 1.2 for all the considered phase space selections. Thus, we 
find that simulations reproduce the shape of the ratio of pf to pj distributions better than the 
individual pf or pj distributions in all selections considered. These four selections mimic phase 
space regions of interest for searches of physics beyond the standard model. We emphasize 
that the agreement is similar for different jet multiplicities and Hj ranges because Z + jets and 
7 -|- jets events have been generated with the same level of accuracy for up to four partons in 
the final-state ME. In the comparison, we considered both processes at either LO or at NLO. It is 
clear from the differences observed between the NLO and LO+PS estimations in each process, 
the conclusions may not be true if the samples are generated with different orders of accuracies 
of the matrix element calculation. 

Our results show that properties of the Z —> pv process can be predicted using the measured 
7 -h jets final state and the simulated ratio between Z ^ pi7 -|- jets and 7 -|- jets. However, this 
simulated ratio must be corrected with the measured ratio of leptonic Z -|- jets and 7 -|- jets. 
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